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Introduction 


Various  military  rockets  are  launched  from  tubes.  The  designer  of 
a  tube-launched  rocket  system  must  consider  the  possibility  of  unbalanced 
forces  acting  on  the  rocket  that  are  caused  by  flow  in  the  annular  gap 
between  the  rocket  and  the  launcher  wall.  As  the  underexpanded  nozzle 
flow  accelerates  upon  leaving  the  nozzle,  it  entrains  air  from  the  annular 
gap  forming  a  shear  layer  at  the  plume  boundary.  The  exhaust  plume  im¬ 
pinges  on  the  wall  a  short  distance  downstream  of  the  nozzle  exit-plane, 
creating  an  impingement  shock-wave.  When  the  impingement  shock-wave  is 
very  weak,  a  fraction  of  the  entrained  air  has  been  given  sufficient  momen¬ 
tum  so  that  it  passes  through  the  impingement  shock  and  the  system  acts  as 
an  ejector.  Using  the  atmosphere  at  the  upstream  end  of  the  launcher  as 
its  source,  a  low  velocity  flow  develops  in  the  annular  regions  to  supply 
the  mass  flow-rate  of  air  that  is  carried  through  the  shock  wave.  A  frac¬ 
tion  of  the  air  entrained  by  the  exhaust  plume  does  not  gain  sufficient 
momentum  to  pass  through  the  shock  wave  and  is  turned  upstream,  creating 
a  separation  "bubble"  near  the  launcher  wall.  Thus,  as  indicated  in  the 
sketch  of  Fig.  la,  the  discriminating  streamline  can  be  traced  into  the 
annular  region  between  the  rocket  and  the  launcher  wall. 

If  the  stagnation  pressure  is  increased  above  the  values  for  ejector- 
type  flow  or  if  the  nozzle  exit-angle  is  relatively  large,  the  flow  at  the 
plume  boundary  must  turn  through  a  larger  angle  when  it  encounters  the 
wall.  Thus,  although  classified  as  a  weak  solution,  the  pressure  rise 
across  the  shock  wave  increases.  Once  the  impingement  shock-wave  becomes 
sufficiently  strong,  the  resultant  adverse  pressure  gradients  causes  a 
fraction  of  the  exhaust  flow  to  be  turned  upstream  into  the  annular  gap, 
as  shown  in  Fig.  lb.  The  exhaust  that  is  turned  upstream  is  known  as  blow- 
by  flow. 

When  conditions  are  such  that  the  impingement  shock-wave  system  is 
strong,  massive  blow-by  occurs.  The  essential  elements  of  a  flow  model 
for  this  launch-tube  flow  field  have  been  identified  experimentally^.  A 
one-dimensional,  theoretical,  flow-model,  which  correlates  well  with  the 
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experimental  data,  has  been  developed  .  The  pressure  gradients  and 
asymmetries  resulting  from  turning  a  large  fraction  of  the  exhaust  flow 
upstream  into  the  annular  gap  cause  unbalanced  forces  on  the  rocket.  Be¬ 
cause  the  rocket  is  flying  free  in  the  launch  tube  at  this  time,  these 
forces  can  cause  the  rocket  to  deviate  significantly  from  its  intended 
flight  path.  Such  deviations  in  the  rocket's  trajectory  have  been  ob- 
served  in  flight  tests’5  and  in  cold-gas  simulations^  .  Therefore,  al¬ 
though  it  is  necessary  to  understand  launch-tube  flow  fields  with  massive 
blow-by,  they  should  be  avoided  whenever  possible. 

To  predict  the  magnitude  and  the  direction  of  the  flow  in  the 
annular  gap,  one  must  be  able  to  describe  the  exhaust  plume  of  the  rocket 
and  the  viscous/shock-interaction  structure  that  results  when  the  plume 
encounters  the  launcher  wall.  The  strength  of  the  impingement  shock  wave 
and  the  characteristics  of  the  viscous  interaction  at  the  wall  depend  on 
the  structure  of  the  exhaust  plume  and  on  the  geometry  of  the  launch  tube. 
When  an  underexpanded,  supersonic  nozzle  exhausts  into  a  constant-area 
tube,  the  strength  of  the  impingement  shock-wave  depends  on  the  Mach  num¬ 
ber  of  the  inviscid  flow  along  the  inner  edge  of  the  mixing  zone  at  the 
jet  boundary,  on  the  ratio  of  the  specific  heats  (y),  on  the  inclination 
of  the  impinging  flow  relative  to  the  launcher  wall,  and  on  the  velocity 
profile  in  the  mixing  zone. 

For  a  wide  range  of  test  conditions,  the  impingement  shock-wave 
is  relatively  weak  and  the  downstream  flow  remains  supersonic.  For  such 
cases,  there  is  relatively  little  secondary  flow  and  the  pressure  in 
the  annular  region  differs  only  slightly  from  the  atmospheric  value.  As 
a  result,  unbalanced  forces  on  the  missile  are  relatively  small.  The 
Chapman-Korst  model6  can  be  used  to  delineate  the  inviscid  and  the  viscid 
flow-mechanisms  encountered  in  the  plume- impingement  flow-field.  Once  the 
geometry  of  the  corresponding  inviscid  plume-boundary  is  established  by  an 
appropriate  potential  flow  analysis,  a  viscous  mixing  component  is  super¬ 
imposed  on  the  boundary.  Wake  closure  conditions  are  specified  in  the  form 
of  a  modified  concept  of  incomplete  turning  of  the  plume  near  the  point  of 
shear-layer  reattachment^. 

Theoretical  flow-fields  have  been  calculated  for  various  nozzle/ 
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launch-tube  configuration  using  this  flow  model.  The  solutions,  thus 
calculated,  are  compared  with  experimental  data  in  the  present  paper. 

The  principal  objectives  of  the  present  paper  are  (1)  to  present  an 
engineering  model  for  the  impingement  flow-field  which  is  produced  when 
an  underexpanded  supersonic  nozzle  is  exhausted  into  a  constant-area  tube, 
(2)  to  compare  the  computed  flow-field  solutions  with  experimental  re¬ 
sults,  and  (3)  to  use  the  comparisons  to  establish  the  validity  of  the 
flow  model . 


Nomenclature 


Cross-section 
Crocco  number 
Diameter 
Mach  number 
Mass  flow- rate 
Pressure 

Radius,  measured  perpendicular  to  the  centerline  of  the  launch  tube 
Initial  radius  of  curvature  of  the  inviscid  plume  boundary 
Temperature 

Streamwise  velocity  component 

Axial  distance  along  the  centerline  of  the  launch  tube  or  streamwise 
distance  (arc  length)  along  the  inviscid  plume  boundary 

Normal  distance  to  the  inviscid  plume  boundary 

Ratio  of  specific  heats 

Nondimensional  shear-layer  coordinate 

Flow  inclination  relative  to  the  centerline  of  the  launch  tube 
Spreading  parameter,  defined  by  equation  (1) 

Dimensionless  velocity  ratio 
Prandtl -Meyer  angle 

Subscripts 

Denotes  static  properties  evaluated  at  the  upstream  end  of  the 
annular  gap 

Denotes  static  properties  in  the  base  region 
Discriminating  (dividing  or  stagnating)  streamline 
Flow  through  the  nozzle  exit-plane,  i.e.,  the  exhaust  flow 
Free-jet  surface  (inviscid  plume  boundary) 

Denotes  properties  evaluated  at  the  intersection  of  the  inviscid 
plume  boundary  and  the  launcher  wall 

Zero  streamline,  i.e.,  that  dividing  the  exhaust  flow  from  the 
base  air 

Denotes  local  conditions  at  the  lip  of  the  nozzle  exit-plane 
Denotes  properties  evaluated  in  the  nozzle  exit-plane 


ne 


5 


Subscripts  (cont.) 

pit  Pressure  measured  using  a  pitot  probe  collinear  with  the  launcher- 
tube  axis 

pk  Denotes  peak  values 

r  Outer  surface  of  the  rocket 

t  Inner  surface  of  the  launch  tube 

tl  Denotes  properties  evaluated  in  the  stagnation  chamber  of  the 
simulated  rocket 

t2  Denotes  properties  downstream  of  a  normal  shock  wave 
v  Beginning  of  viscous  interaction 

1  Mass  flow-rate  given  by  equation  (6) 


Theoretical  Analysis 

A  component  approach  for  supersonic,  turbulent  flows,  i.e.,  the 
Chapman-Korst  model^,  forms  the  basis  for  the  present  analysis.  It 
includes  the  delineation  of  inviscid  and  viscid  flow-mechanisms  encountered 
in  plume/launcher-wall  impingement  flows  and  their  subsequent  synthesis 
into  an  overall  systems  solution.  The  geometry  of  the  corresponding 
inviscid  plume  boundary  is  established  by  an  appropriate  potential  flow 
analysis.  The  viscous  mixing  component  is  superimposed  on  this  boundary. 
Wake  closure  conditions  are  specified  in  the  form  of  a  modified  concept  of 
the  incomplete  turning  of  the  plume  near  the  point  of  shear  layer  reattach¬ 
ment'7.  A  number  of  conceptual  simplifications  for  the  plume-impingement 
flow  field,  which  will  be  introduced,  are  justifiable  for  the  conditions  of 
i nterest. 

Since  the  launching  of  rockets  from  straight  tubes  will  generally 
impose  geometric  restrictions  to  blow-by,  the  present  analysis  can  be  sim¬ 
plified  by  the  assumption  that  the  mass  flow-rates  which  will  be  rejected 
from  or  entrained  into  the  viscous  jet  mixing  region  should  be  small. 

This,  in  turn,  allows  the  adoption  of  a  linear  velocity  profile  in  the 

g 

mixing  region  .  It  is  assumed  that  the  entire  system  will  be  operating 
under  essentially  isoenergetic  conditions,  that  is,  with  the  stagnation 
temperature  equal  to  that  of  the  exhaust  gas.  This  assumption  is  realistic 
for  most  flows  under  consideration;  the  exception  being  ejector-type  flow 
generated  by  an  actual  rocket  firing. 

Since  the  ratio  of  the  launch-tube  diameter  to  rocket  nozzle-exit 
diameter  is  limited  to  values  £  1.5,  the  contour  of  the  corresponding 
inviscid  plume  boundary  can  be  determined  with  sufficient  accuracy  by  the 
second-order  potential -flow  approximation  for  axially  symmetric,  Sj.;;er- 

Q  IQ 

sonic  flow  near  the  center  of  an  expansion  ’.  The  resulting  circu'at- 
arc  configuration  is  convenient  for  establishing  a  simple  inviscid-flow- 
component  geometry.  On  the  other  hand,  it  is  necessary  to  treat  the 
nozzle  flow  in  detail  in  order  to  define  accurately  the  conditions  at,  and 
near,  the  lip  of  the  exit  cross-section.  This  requires  the  analyses  of  the 
nozzle  flow  by  the  method  of  characteristics  (including  the  effects  of 
characteristics  coalescence)  as  it  develops  from  the  transonic  flow  region 
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near  the  nozzle  throat^. 

The  exit-plane  conditions  near  the  nozzle  lip,  will  be  determined 
by  the  local  wall-slope,  6^,  the  local  wall-curvature,  the  local  Mach  num¬ 
ber,  M.  ,  and  the  local  flow  acceleration.  With  these  initial  conditions 
specified,  selection  of  the  jet  surface  Mach  number,  Mp,  will  yield  a 
unique  solution  for  the  circular-arc  approximation  of  the  contour  of  the 

inviscid  plume  boundary  in  terms  of  the  initial  streamline  angle,  6r,  and 

12  ^ 

the  radius  of  curvature, 

Thus,  for  a  given  nozzle  configuration  and  for  any  jet  surface  Mach 

number,  Mp,  the  inviscid  flow  field  will  be  uniquely  determined,  see  Fig. 

2,  with  0.  denoting  the  inviscid  impingement  angle  at  the  cylindrical 

launcher  wall,  i.e.,  at  radius  r^.  One  now  has  to  analyze  the  viscous  flow 

mechanisms  which  will  provide  the  closure  conditions  for  determining  the 

actual  operating  condition  of  the  system. 

As  mentioned  above,  the  jet  mixing  characteristics  of  the  free  shear- 

layer  at  the  plume  boundary  can  be  approximated  by  selecting  isoenergetic 

conditions  and  a  linear  velocity  profile*.  The  spread  parameter  o  shall 
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here  be  related  to  the  rate  of  change  in  slope  of  the  velocity  profile  ’ 
in  the  form: 


o 


1 


—  [(— )  1 

Dx  LV9y;maxJ 


(1) 


For  a  linear  velocity  profile  with  the  conditions  at  the  "inviscid"  plume 
boundary  given  by  Up  and  ip,  one  finds: 


where  n 


(2) 


The  density  ratio  in  the  iso-energetic  mixing  region  is: 


★ 

Use  of  a  self-similar,  linear  mixing  profile  does  not  exclude  the  study  of 
Reynolds  number  effects.  The  effects  of  the  boundary  layer  along  the 
nozzle  wall  can  be  considered,  subsequently  either  in  form  of  an  equiva- 
7  13 

lent  bleed  ’  or  by  analyzing  the  development  of  a  dissipative  shear 
layer  in  a  rotational  but  essentially  non-dissipative  flow  field 
envisioned  by  the  flow  expansions  near  the  nozzle  lip. 
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PF 


1  -  C 
1  -  o2C 


where  use  of  the  Crocco  number. 


CF  =  +  MF2 


(3) 


(4) 


conveniently  eliminates  the  influence  of  the  specific  heat  ratio  on  con¬ 
stant-pressure  mixing  profiles. 

Application  of  the  momentum  equation  to  a  two-dimensional  mixing 
region  allows  one  to  locate  the  profile  with  respect  to  the  corresponding 
inviscid  jet  boundary  in  the  form: 


nM  =  °  x~  = 


1  - 

[i  -  c2 

_2 

V 

l  c 

-1  +  2C  In 


1  +  C 


(5) 


One  should  note  that  the  application  to  an  axi symmetric ,  curved  plume  boun 
dary  requires  that  the  length  be  interpreted  as  an  arc  length  and  that 
the  use  of  o  values  originally  defined  for  planar  mixing  zones  adds  to 
the  inherent  uncertainty  in  making  an  assessment  for  the  effects  of  com¬ 
pressibility  on  the  spread  parameter.  For  the  present  purpose,  where  the 

Mach  number  range  for  MF  will  seldom  exceed  3.5,  one  may  use  the  simple 

16  ' 
relationship  : 


o  =  12  +  2.56  MF 

The  mass  rate  of  flow  per  unit  depth  within  the  mixing  zone  up  to  the 


streaml ine  y: 


,yU) 


,:,1  = 


pu  dy 


y+-o 


can  now  be  expressed  in  closed  form  by: 


m1  1 1  ( 4>  »Cp) 


x  u 


F  MF 


f  I  ■  IMjjjlliiiig 


a 


(6) 
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where  I^(<J>,Cp)  =  /r T" 


cf2-i 


In  (1  -  /cp2)  ; 


the  velocity  ratio  <p  -  —  defines  a  streamline  having  the  location 


y  =  —  (nM  -  — -■)  measured  perpendicular  to  and  outside  of  the  correspond¬ 
ing  inviscid  plume  boundary.  The  velocity  along  the  zero  streamline,  which 
separates  the  nozzle  flow  from  the  air  entrained  from  the  base  region  is 
given  by 

♦j  ■  CF2)  Exp  ^ ,n  (r^|  -2]  <7> 

The  recompression  of  the  plume  as  it  impinges  on  the  cylindrical 
wall  leads  to  a  flow  configuration  as  shown  in  Fig.  2b.  If  the  stagnating 
(discriminating)  streamline  has  a  velocity  level,  <|)d,  different  from  that 
of  the  zero  streamline,  <j>.,  a  secondary  flow  will  be  established  as  re- 
quired  for  satisfying  the  conservation  of  mass  in  the  near  wake  region. 

The  secondary  flow  rate  will  be  positive,  i.e.,  ejector-type  flow,  if 
<{>j  >  and  negative,  i.e.,  blow-by,  if  <pj  <  <f>^.  The  secondary  flow  rate 
may  be  approximated  as 

m  -  ^nr^  —  UpPp[1 1  ( <T>  j  ,Cp )  "Ii(<t>j,Cp)J  (8) 

The  condition  <Jk  =  corresponds  to  an  operation  without  any  secondary 
flow  which  will  be  satisfied  for  near  wake  vented  to  the  atmosphere  if 

p  _  patm 

^  „  '  pti 

F 

or  the  case  where  the  annular  gap  is  isolated  by  a  tight  fitting  sabot  or 
by  an  obturator  ring. 

As  the  discriminating  streamline,  whose  velocity  ratio  is  <j>.,  stagnates 

u  6  7 

at  the  wall,  the  process  of  recompression  will  be  practically  isentropic  ’  . 
Hence,  for  a  local  wall  pressure  corresponding  to  an  inviscid  value  of 
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the  Crocco  number,  Cd  =  Cp<j>d, 


Since  the  reattachment  point  of  this  streamline  falls  into  a  region 
where  the  adjacent  inviscid  plume  flow-field  has  not  yet  fully  aligned 
with  the  solid  wall,  the  concept  of  an  incomplete  turning  angle  will  be 
utilized. 

For  a  linear  velocity  profile,  the  correlation  proposed  by  Page,  et 
al.^  for  cases  involving  mass  bleed  reduces  to  a  remarkably  simple  ex- 
PreSSi0n:  U)(Cr)  -  ia(Cj 

A  ■  ^T^cfy  --  °-59  <10) 

which  is  valid  over  a  range  of  inviscid  boundary  Mach  number  from  1  to  5, 
i.e.  1  <_  Mp  _<  5.  In  equation  (10),  w(C)  is  the  Prandtl-Meyer  function 
(which  depends  on  the  Crocco  number)  and  is  the  Crocco  number  for  the 

inviscid  flow  after  the  impingement  of  the  inviscid  jet  and  its  realignment 
with  the  launcher  wall. 

Equations  1  through  10  can  be  solved  to  establish  the  theoretical 
operating  conditions  using  information  about  the  base  pressure,  which  is 
approximately  equal  to  the  atmospheric  value  when  the  annular  gap  is  vented, 
and  the  peak  reattachment  pressure,  which  is  located  a  distance: 

xi  =  R^sin  eL  -  sin  9^  (11) 

from  the  nozzle  exit-plane.  For  cases  where  the  secondary  flow  rates  are 
large,  equations  1  through  10  must  be  supplemented  by  a  relationship  be¬ 
tween  the  mass  flow-rate  and  the  pressure  drop  in  the  annular  gap.  This 
additional  relation  for  the  large  flow-rates  in  the  annular  gap  is  dis¬ 
cussed  in  Ref.  2. 


The  inviscid  pressure  rise  will  start  near  the  location,  see  Fig.  2b, 
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where  the  implicit  equation  for  6 
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has  to  be  solved  by  iteration. 

After  plume  impingement  and  alignment  of  the  jet  surface  with  the 
wall  by  reversible  compression,  the  pressure  distribution  for  the  inviscid 
flow  configuration  can  again  be  treated  with  reasonable  approximation  by 
the  formal  application  of  the  concept  of  centered  expansions.  With  the 
new  origin  of  the  centered  expansion  of  the  impingement  point  for  the 
inviscid  plume,  an  ordinary  differential  equation  is  obtained  for  M*  along 
the  wall  downstream  of  the  impingement  and  the  resulting  pressure  distri¬ 
bution  can  be  determined.  The  mathematical  approximations  inherent  in  the 
method  of  centered  expansions  restricts  the  useful  range  of  the  calculated 
wal 1 -pressure  distribution  to  approximately  1.0  r^ downstream  of  the 
impingement. 

It  is  of  interest  to  explore  the  influence  of  Reynolds  number  to¬ 
gether  with  an  assessment  of  the  nozzle-wall  boundary-layer  development. 

The  boundary  layer  in  axisymmetric,  convergent/divergent  nozzles  can  be 
18  19  20 

calculated  ’  ’  at  the  nozzle  lip.  Then,  for  a  given  Mach  number,  Mp, 

one  can  calculate  the  development  of  the  dissipative  shear  layer  follow¬ 
ing  the  expansion  from  the  nozzle  exit-plane.  After  the  transition  of  the 
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thin  laminar  sublayer  ,  the  turbulent  mixing  layer  grows  until  its  inner 
edge  reaches  and  then  crosses  into  the  "free-stream"  within  the  rotational 
shear  layer  produced  by  the  expanded  nozzle  boundary  layer.  The  develop¬ 
ing  mixing  layer  proceeds  toward,  but  usually  does  not  reach,  self-similar 
mixing  profiles.  Calculations  for  the  present  test  conditions  indicate  that 
the  velocity  along  the  dividing  streamline  is  within  five  percent  of  the 

similarity  value.  Thus,  the  test  simulations  of  the  current  program  provide 

a  reasonable  simulation  for  the  prototype  reattachment  pressure  rise. 

Of  greater  concern  is  the  behavior  of  the  reattachinq  flow  due  to 
the  strong  interaction  with  its  boundary  layer.  Implied  with  the  concept 
of  incomplete  turning  of  the  adjacent  inviscid  stream  is  the  continuation 
of  the  pressure  rise  beyond  the  stagnation  point  within  the  dissipative  layer. 
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This  will  lead  to  a  thickening  of  the  attached  boundary  layer  and  may,  by 
itself,  cause  a  pressure  overshoot  above  the  level  consistent  with  the  full 
turning  of  the  external  flow. 

On  the  other  hand,  there  will  also  be  a  communication  within  the 
boundary  layer  of  the  rapid  expansion  following  the  plume  impingement. 

These  mechanisms  are  expected  to  modify  the  pressure  distribution  across 
the  reattachment  region,  but  in  opposing  fashion.  As  a  result,  pressure 
peaks  determined  for  the  corresponding  inviscid  jet  boundary  will  generally 
show  reasonable  agreement  with  experimental  data  even  though  there  may  be 
noticeable,  albeit  explainable,  differences  between  theoretically  deter¬ 
mined  and  measured  pressure  distributions. 


Experimental  Program 


The  experimental  results  presented  in  this  paper  were  obtained  at 
the  Rocket  Exhaust  Effects  Facility  located  at  the  Experimental  Aerody¬ 
namics  Laboratory  (EAL)  of  the  University  of  Texas  at  Austin.  Simulated 
rocket  exhaust  plumes  were  obtained  by  accelerating  unheated,  compressed 
air  (the  test  gas)  through  a  convergent/divergent  nozzle  (the  simulated 
rocket).  The  nozzle  flow  was  exhausted  either  into  quiescent  air  or  into 
a  constant-area  launch  tube.  A  sketch  of  a  typical  nozzle/launch-tube  con¬ 
figuration  is  presented  in  Fig.  3.  The  simulated  launch  tubes  were  mounted 
on  a  moveable  table,  so  that  the  location  of  the  launch  tube  relative  to 
the  exit  plane  of  the  simulated  rocket  nozzle  could  be  varied.  For  all  of 
the  runs  in  the  present  test  program,  the  axis  of  the  launch  tube  was 
col  linear  with  the  axis  of  the  simulated  rocket  nozzle. 

The  test  schedule  included  three  different  nozzle  geometries  and 
three  launch  tubes.  Since  the  objective  of  the  present  paper  is  to  use  the 
experimental  data  to  validate  the  theoretical  flow  model  and  to  establish 
the  limits  of  its  validity,  only  selected  data  are  presented  in  this  paper. 
For  a  more  complete  discussion  of  the  experimental  programs,  the  reader  is 
referred  to  refs.  22  and  23. 


Simulated  Rockets 

Simulated  rocket  exhaust  flows  were  produced  by  accelerating  unheated 
air  through  convergent/divergent  nozzles.  Data  are  presented  for  two  of  the 
nozzle  configurations, 

1)  a  20°  conical  convergent  divergent  nozzle,  designated  as  C3,  and 

2)  a  10°  conical  contoured  convergent  divergent  nozzle,  designated 
as  C5. 

The  geometries  of  the  nozzle  configurations  studied  in  the  current 
program  are  presented  in  the  equations  below. 

The  20°  Conical  Nozzle,  C3: 

For  -1.308  ^  x  <_  -1.031: 

(x  +  1.308)2  +  (r  -  1.626)2  *  (0.813)2 
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For  -1 .031  £  x  £  0.0: 

r  =  0.759  +  0.364  (x  +  1.308) 

The  10°  Contoured  Nozzle,  C5: 

For  -2.908  £  x  £  -1.895: 

(x  +  2 . 098 ) 2  +  (r  -  1.626)2  =  (0.813)2 
For  -1 .895  £  x  £  0.0: 
r2  =  1.530  +  0.434x 

The  equations  give  the  dimensions  in  centimeters.  A  sketch  of  the  C5  nozzle 

is  presented  in  Fig.  4.  The  throat  radius,  the  nozzle  exit-plane  radius, 

and  the  external  radius  of  the  rocket  were  the  same  for  all  nozzles. 

The  cross-sectional  area  for  all  nozzle  exit-planes  was  A  =  2.316A*. 

r  ne 

If  the  acceleration  of  the  flow  through  the  nozzles  of  this  area  ratio  were 
isentropic,  the  Mach  number  of  the  nozzle  exit  plane  would  be  2.36  (ref.  24). 


Simulated  Launch  Tubes 

The  constant  area  launch  tubes,  designated  as  LI,  L2,  and  L3,  were 
used  for  the  tests  (Fig.  3).  The  inside  radius  (rfc)  for  the  smooth  bored 
tubes  was  1.20  rng  (1.484  cm)  for  the  LI  configuration,  1.275  rne  (1.522  cm) 
for  the  L2  configuration,  and  1.50  rng  (1.855  cm)  for  the  L3  configuration. 
Static  pressure  orifices  (0.0635  cm  in  diameter)  were  located  axially  along 
the  tube  wall  at  a  spacing  of  0.25  rne  (0.309  cm).  Additional  static 
pressure  orifices  spaced  0.50  rne  (0.618  cm)  apart  were  located  further 
downstream. 

For  each  test  configuration,  i.e.,  nozzle/launch-tube  configuration, 
two  pitot-static  probes  were  located  midway  between  the  simulated  rocket 
and  simulated  launcher  wall  and  4  1 28  cm  from  the  upstream  end  of  the  launch 
tube.  Static  pressure  orifices  were  located  in  the  launch-tube  wall  at  the 
face  of  the  pitot  probes.  These  pitot-pressure/static-pressure  data  could 
be  used  to  "measure"  the  mass  flow-rate  in  the  annular  gap.  The  pitot  probe 
was  designed  so  that  it  could  face  either  upstream  or  downstream  in  order 
to  "measure"  the  mass  flow-rate  in  either  direction,  i.e.,  ejector  flow  or 
blow-by  flow. 
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The  Mass  Flow-Rate 

The  mass  flow-rates  were  calculated  using  the  pressure  data  from  the 
pitot-static  probes  located  at  the  upstream  end  of  the  annular  gap.  If  one 
assumes  one-dimensional,  isentropic  flow  in  the  annular  gap,  and  perfect 
gap  properties,  then  the  dimensionless  mass  flow-rate  is  given  by: 


For  those  runs  where  blow-by  occurred,  p^„  was  assumed  to  be  the 

atmospheric  value  and  the  measured  pitot  pressures  were  used  to  calculate 

M  .  For  the  ejector  flow  calculations,  the  pitot  pressures  were  assumed 
a9 

to  be  the  atmospheric  value  while  the  measured  value  of  pflq  was  used  in 
the  calculation  of  the  Mach  number  in  the  annular  gap. 


Obturator  Rings 

For  some  of  the  runs,  obturator  rings  were  used  to  seal  the  base 
region  so  that  there  could  be  no  flow  in  the  annular  gap.  The  obturator 
ring  fit  snugly  into  the  annular  gap  at  an  x  of  approximately  -1.770  cm 
i.e.,  upstream  of  the  nozzle  exit  plane. 


Discussion  of  Results 


Solutions  of  the  theoretical  flow-fields  have  been  computed  from 
the  transonic  region  of  the  nozzle  throat,  through  the  nozzle,  and  through 
the  impingement  of  the  jet-exhaust  plume  on  the  launcher-wall.  Computed 
solutions  have  been  obtained  both  for  the  C3  and  for  the  C5  nozzles. 

Thus,  the  theoretical  solutions  of  the  rocket  exhaust  flow-field  represent 
the  effects  of  the  actual  nozzle  geometries  on  the  overall  flow-field.  The 
theoretical  solutions  will  be  compared  with  the  experimental  data.  In 
some  cases,  two  theoretical  solutions  representing  varying  degrees  of 
rigor  were  generated.  In  such  cases,  the  theoretical  solutions  will  be 
compared  with  each  other. 


Nozzle  Flow  Field 

The  static  pressures  were  measured  along  the  wall  of  the  divergent 
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section  of  the  C5  nozzle  for  stagnation  pressures  from  0.7  x  10  N/m  to 

f.  ?  fi  ? 

8.8  x  10  N/m  .  For  stagnation  pressures  of  3.4  x  10  N/m,  or  greater, 
the  values  of  the  ratio  P/Ptl  were  a  function  of  position  only.  The  non- 
dimensional  ized  static-wal 1 -pressure  measurements  are  compared  in  Fig. 

5  with  the  theoretical  values  computed  using  the  method  of  characteristics. 
The  agreement  between  the  theoretical  and  the  experimental  distributions 
is  excellent. 

The  radial  pitot-pressure  distributions  which  result  when  unheated 
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air  is  accelerated  from  a  reservoir  where  p^  is  8.715  x  10  N/m  through 
the  C3  nozzle  and  exhausts  into  quiescent  air  are  presented  in  Fiq.  6. 

The  theoretical  values  of  p^/p^  are  compared  with  experimental  pitot 
pressure  for  the  nozzle  exit-plane  (x  =  0.0  rne)  and  for  the  plane  at 
x  =  0.5  rne.  Theoretical  values  are  presented  only  for  the  core  reqion. 

The  theoretical  values  of  the  pitot  pressures  were  calculated  using  the 
local  Mach  number,  which  is  computed  using  an  isentropic  method-of- 
characteristics  solution,  to  define  the  stagnation  pressure  ratio  across 
a  normal  shock  wave  ( pt2 ^ fJt  1  ^ the  ^oca^  fl°w  inclination 
was  not  included  in  the  calculation  procedure.  The  effect  on  the 
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theoretical  stagnation  ratios  presented  in  Fig.  6  should  be  small,  since, 
as  Pope  and  Goin  note,  a  pitot  probe  "will  be  insensitive  to  angle  of 
attack  up  to  10  deg  for  an  orifice  diameter  of  only  10  percent  of  the 
outside  diameter  and  up  to  15  deg  for  one  98  percent  of  the  outside 
diameter."  The  theoretical  total -pressure  distributions  in  the  nozzle 
exit-plane  are  in  very  good  with  the  experimentally-determined  distri¬ 
butions.  Referring  to  the  tabulated  values  of  Ref.  24,  the  total  pres¬ 
sures  indicate  that  the  Mach  number  varies  from  2.13  near  the  axis  to 
2.48  near  the  wall  of  the  nozzle.  Recall  that  if  one  assumes  that  the 
flow  undergoes  a  one-dimensional  acceleration  from  the  sonic  throat  to 
the  stream  tube  of  the  exit  plane,  the  Mach  number  would  be  2.36.  Thus, 
although  the  C3  is  a  conical  nozzle,  the  flow  in  the  nozzle  exit-plane  is 
not  conical. 

For  x  =  0.5  rne,  the  pitot-pressure  measured  in  the  internal  region 

are  in  reasonable  agreement  with  the  theoretical  values.  The  theoretical 

model  even  predicts  the  variation  which  occurs  near  r  =  0.5  r  .  For 

ne 

r  >  0.965  rng,  the  measured  pitot-pressures  decrease  rapidly  with  increas¬ 
ing  r  until  near  the  plume  boundary.  These  data  illustrate  that  since  the 
Mach  number  is  constant  along  an  "inviscid  plume  boundary",  the  local 
Mach  number  is  a  maximum  at  a  point  between  the  edge  of  the  core  region 
and  the  plume  boundary. 

O 

When  the  stagnation  pressure  is  8.715  x  10  N/m  ,  the  static  pres¬ 
sure  in  the  nozzle  exit-plane  is  roughly  six  times  the  atmospheric  value. 
Thus,  as  the  exhaust  gases  leave  the  nozzle,  they  accelerate  rapidly 
until  the  pressure  along  the  plume  boundary  is  equal  to  the  static  pres¬ 
sure  of  the  surrounding  gas.  The  plume  boundaries  for  the  C3  (20°  coni¬ 
cal)  and  for  the  C5  (10°  contoured)  nozzles,  as  calculated  using  the 
method-of-characteristics  are  presented  in  Fig.  7.  The  plume  boundaries 
computed  using  the  method-of-characteristics  are  in  excellent  agreement 
with  those  observed  in  the  schlieren  photographs  of  the  nozzles  exhausting 
into  quiescent  air.  This  is  true  for  both  of  these  nozzles,  which  are  of 
considerably  different  geometries. 
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The  correlation  between  the  theoretical  and  the  experimental ly- 
determined  flow-fields  indicates  that  the  nozzle-exhaust  flow-field  can 
be  well  modeled  analytically.  Thus,  the  method-of-characteristics  for 
axisymmetric  flow  provides  an  acceptable  tool  even  when  no  direct  pro¬ 
visions  are  made  for  the  occurence  of  internal  shocks.  As  will  be  seen, 
the  plume  impingement  and  the  initial  wal 1 -pressure  rise  will  be  described 
with  sufficient  accuracy  by  the  method-of-characteristics  followed  by 
plane  shock  relations. 

The  internal  radii  of  the  launch  tubes  used  in  the  present  test 
program  are  from  1.2  rpe  to  1.5  rng,  depending  on  the  launch-tube  con¬ 
figuration.  Since  the  developing  shear  layer  is  relatively  short,  the 
ensuing  viscous/shock-interaction  structure  at  the  wall  is  very  dependent 
on  the  expansion  process  as  the  flow  leaves  the  nozzle. 

For  the  relatively  small  values  of  rt/rne  of  interest  in  tube- 
launched-rocket  applications,  it  is  necessary  to  model  the  plume  boundary 
only  for  a  short  distance.  Thus,  the  calculation  of  the  inviscid  plume 
boundary  by  a  technique  which  provides  quick,  but  reasonable  contours  be¬ 
comes  attractive.  The  theoretical  plume  boundaries  have  been  calculated 

9  10 

using  the  method-of-centered-expansions  ’  .  These  approximate  plume 

boundaries  are  compared  with  those  calculated  using  the  more  rigorous 
method-of-characteristics  and  with  the  data  in  Fig.  7.  Note  that  the 
plume  boundary  calculated  using  the  method-of-centered-expansions  is  very 
close  to  the  actual  boundary  for  both  nozzles.  Flowever,  at  a  given  x- 
location,  the  local  flow  angle  at  the  plume  boundary  calculated  using  the 
method-of-centered-expansions  is  slightly  less  than  the  actual  flow  angle. 
Thus,  the  inviscid  plume  boundaries  calculated  using  the  method-of- 
centered-expansions  provide  a  close  approximation  of  the  actual  plume 
boundaries,  making  the  method  attractive  for  use  in  generating  quick  calcu¬ 
lations  of  the  plume  impingement  flow. 
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Launch-Tube  Flow  Fields 

When  a  highly  underexpanded  supersonic  nozzle  exhausts  into  a  con¬ 
stant-area  tube,  the  viscous/shock-wave  interaction  due  to  the  impinging 
flow  is  subject  to  a  boundary  condition  in  the  base  region.  If  the  annular 
gap  between  the  rocket  and  the  launch-tube  wall  is  open,  the  atmospheric 
pressure  at  the  forward  end  of  the  vented  annular  gap  forms  the  boundary 
condition.  The  direction  and  the  magnitude  of  the  secondary  flow  which  is 
usually  generated  by  the  impingement  interaction  is  controlled  by  this 
boundary  condition.  If  the  mass  flow-rate  of  the  secondary  flow  in  the 
annular  gap  is  relatively  low,  the  pressure  will  vary  only  slightly  along 
the  annular  gap.  Thus,  the  static  pressure  in  the  base  region  will  be 
approximately  equal  to  the  atmospheric  value.  However,  if  the  annular 
gap  is  sealed,  the  pressure  in  the  base  region  will  assume  a  value  such 
that  the  viscous/shock-wave  interaction  at  impingement  creates  a  "constant- 
mass"  recirculating  flow  in  the  base  region. 

Pressure  distributions  for  the  launch-tube  flow-fields  with  the  annular 
gap  sealed  by  an  obturator  ring.  -  Over  a  wide  range  of  stagnation  pres¬ 
sure,  the  ratio  p^/p^  is  essentially  independent  of  the  stagnation  pres¬ 
sure  when  the  annular  gap  is  sealed.  That  the  pressure  ratio  p^/p^ 

should  be  constant  for  a  given  conf iguration  for  a  fully  developed  super- 
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sonic  flow  is  well  documented,  e.g.,  refs.  23,  26,  and  27,  Korst  et  al 
have  shown  that  the  equilibrium  (or  steady-state)  value  of  the  base 
pressure  is  governed  by  the  viscous  shear-layer  at  the  boundary  of  the 
nozzle  exhaust  plume  that  bounds  the  "dead-air"  base-flow  region  and  by 
the  physical  geometric  constraints.  The  theoretical  values  of  the  base 

pressure  are  presented  for  the  C3  and  the  C5  nozzles  in  Fig.  8  for  all 

three  launch  tubes  for  stagnation  pressures  from  2.5  x  10°  N/rn  to 

8.8  x  1G6  N/m2. 

The  experimental  values  of  the  static  wall-pressure  vary  only 
slightly  with  x  in  the  region  between  the  obturator  ring  (small  negative 
x-coordinate)  and  the  intersection  of  the  exhaust  plume  with  the  launcher 
wall,  reaching  a  minimum  just  before  the  impingement-induced  pressure  rise. 
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This  is  illustrated  by  the  data  presented  in  Figs.  9  and  10.  For  the 
purposes  of  the  present  paper,  the  static  wall -pressure  sensed  by  the  ori¬ 
fice  located  at  x  =  0.0  rn0  (i.e.,  in  the  plane  of  the  nozzle  exit)  is 
designated  as  the  base  pressure.  The  experimental  values  of  pb  are  pre¬ 
sented  in  Fig.  8.  Since  they  are  a  linear  function  of  p  ^ ,  the  ratio 
Pb/Pti  is  constant  for  a  given  configuration.  Furthermore,  they  compare 
very  well  with  the  theoretical  values  for  both  nozzles. 

The  static-pressure  distributions  in  the  vicinity  of  the  plume 
impingement  have  been  computed  for  flows  with  the  annular  gap  sealed  by 
an  obturator  ring.  The  pressure  distributions  for  unheated  air  accelera- 
ting  from  a  reservoir  where  p  ^  is  8.715  x  10  N/m  through  the  C3  nozzle 
into  the  constant-area  launch-tubes  are  presented  in  Fig.  9.  Similar 
solutions  for  the  C5  nozzle  are  presented  in  Fig. 10.  The  theoretical 
solutions  are  compared  with  experimental  data  for  these  nozzle/launch- 
tube  configurations.  Note,  that  as  the  tube  radius  is  increased  from 
1.20  rne  (the  LI  tube)  to  1.50  rng  (the  L3  tube),  the  base  pressure  de¬ 
creases  and  the  location  of  the  peak  pressure  moves  further  downstream 
of  the  nozzle  exit-plane.  However,  as  the  base  pressure  decreases,  the 
angle  through  which  the  plume  boundary  expands  at  the  nozzle  exit-plane 
increases.  As  a  result,  the  pressure  ratio  across  the  impingement  shock- 
wave  (i.e.,  the  peak  static-pressure  divided  by  the  base  value)  increases 
as  r^  increases. 

For  both  nozzles,  the  theoretical  value  of  the  peak  pressure  in  the 
impingement  region  is  significantly  greater  than  the  corresponding  experi¬ 
mental  value  for  the  LI  launch-tube.  Because  of  the  small  difference 
between  r^  and  rn0,  the  free-shear  layer  at  the  plume  boundary  is  relatively 
short  and  the  flow  deceleration  associated  with  the  impingement  process 
significantly  modifies  the  vi scous/shock-wave  interaction.  This  modifi¬ 
cation  is  not  reflected  in  the  theoretical  flow  model,  and  thus,  the  peak 
pressure  is  overestimated  by  the  theory.  However,  the  agreement  between 
the  theoretical  and  the  experimental  values  are  much  better  for  the  L3 
launch-tube.  The  improvement  in  the  correlation  for  the  larger-radius 
launch-tube  is  attributed  to  compensating  effects  of  approximations  in  the 
theoretical  flow  model. 
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Pressure  distributions  for  the  launch-tube  flow-fields  with  a  vented 
annular  gap.  -  As  noted  earlier,  if  the  mass  flow-rate  of  the  secondary 
flow  in  the  annular  gap  is  relatively  low,  the  pressure  will  vary  only 
slightly  along  the  annular  gap.  Since  the  static  pressure  in  the  nozzle 
exit-plane  (p  )  is  directly  proportional  to  the  stagnation  pressure, 
while  the  base  pressure  (pb)  remains  roughly  constant,  the  jet  pressure 
ratio  (Png/P^)  varies  with  ptl .  Thus,  the  viscous/shock-wave  interaction 
depends  on  the  stagnation  pressures  as  well  as  the  nozzle/launch-tube 
configuration. 

The  effect  of  the  launch-tube  radius  is  indicated  by  the  pressure 
distributions  presented  in  Figs.  11  and  12  for  the  largest  value  of  p^ 
tested.  Since  the  base  pressure  is  approximately  constant,  the  shape  of 
the  plume  boundary  would  be  independent  of  the  launch-tube  radius,  if  the 
flow  were  inviscid.  If  the  shape  of  the  plume  boundary  were  indeed  inde¬ 
pendent  of  the  launch-tube  radius,  as  r^  decreased,  the  intersection  of 
the  plume  boundary  and  the  wall  would  be  nearer  the  nozzle  exit.  Thus, 
the  impingement  shock-wave  would  move  upstream  and  increase  in  strength 
as  rt  decreases. 

These  trends  are  substantiated  by  both  the  theoretical  and  the 
experimental  pressure  distributions  for  the  nozzle/launch-tube  configura¬ 
tions  using  the  C5  nozzle,  which  are  presented  in  Fig.  12.  As  was  the 
case  for  the  flow  fields  with  the  obturator  ring,  the  correlation  between 
theory  and  experiment  is  better  for  the  larger  radius  launch-tube,  the 
L3.  As  before,  the  theoretical  model  of  the  shear  layer  does  not  reflect 
the  interaction  between  the  approaching  free  shear-layer  and  the  pressure 
rise  of  the  impingement  shock-wave.  However,  as  was  shown  in  Fig.  7,  the 
method-of-centered-expansions,  which  is  used  to  provide  the  approximate 
inviscid  plume  boundary,  yields  plumes  for  which  the  local  flow  angle  at 
the  plume  boundary  is  slightly  less  than  the  actual  flow  angle.  Because 
these  two  flow-model  approximations  tend  to  compensate  each  other,  the 
theoretical  values  of  the  peak  pressure  are  in  very  good  agreement  with 
the  data  for  these  flow/geometry  conditions. 
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Although  experimental  pressure  distributions  are  presented  for  all 
three  C3/launch-tube  configurations,  only  for  the  C3/L3  is  the  theoretical 
distribution  presented.  Theoretical  solutions  were  not  possible  for  the 
C3/L1  or  for  the  C3/L2  configurations  with  the  present  flow  model,  since 
the  local  flow  inclination  at  the  plume  boundary  intersection  with  the 
launcher  wall  exceeds  the  value  which  can  be  turned  through  a  weak  shock- 
wave.  This  contrasts  the  flow  fields  for  the  sealed  annular  gap.  Recall 
that  theoretical  pressure  distributions  were  presented  in  Fig.  9  for  all 
three  launch  tubes  when  the  annular  gap  was  sealed  by  an  obturator  ring. 

As  shown  in  Fig.  8a,  the  base  pressure  is  considerably  greater  than  the 
atmospheric  value,  when  the  annular  gap  is  sealed,  and  the  expanding 
flow  turns  through  a  smaller  angle  as  it  leaves  the  nozzle.  Therefore, 
with  an  obturator  ring  in  place,  the  local  flow  inclination  of  the  imping¬ 
ing  plume  is  such  that  the  impingement  shockwave  is  weak. 

Although  the  theoretical  solutions  indicate  the  weak  impingement 
shock-waves  are  not  possible  for  the  C3/L1  or  the  C3/L2  configurations, 
the  maximum  measured  values  of  the  pressure  are  not  large  enough  to 
indicate  the  presence  of  a  strong  shock-wave.  Thus,  these  pressure  distri¬ 
butions  indicate  a  relatively  strong  interaction  between  the  impingement 
shock-wave  and  the  free-shear  at  the  plume  boundary.  For  the  two  smaller 
radii  tubes,  the  streamwise  pressure-gradient  along  the  launcher  wall  is 
not  constant  in  the  impingement  region.  As  was  discussed  in  ref.  23,  the 
relatively  strong  viscous/shock-wave  interaction  which  occurs  for  this 
nozzle/launch-tube  configuration  at  high  stagnation  pressures  produces 
significant  blow-by  flow. 

Secondary  flows.  -  As  the  rocket  exhaust  entrains  base-region  air,  a 
free-shear  layer  develops  along  the  plume  boundary,  accelerating  air 
from  the  base-region  and  decelerating  exhaust  gas.  When  this  free-shear 
layer  encounters  the  adverse  pressure  gradient  associated  with  the  impinge¬ 
ment  shock-wave,  some  of  the  flow  passes  through  the  shock  wave,  some  of 
it  is  turned  back.  As  discussed  in  Fig.  1,  if  the  stagnation  pressure 
is  relatively  low,  the  system  acts  as  an  ejector  with  some  of  the  entrained 
air  passing  through  the  impingement  shock-wave.  As  the  stagnation  pressure 
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is  increased,  the  ratio  Pne/Pb  increases  and  the  flow  expands  through  a 
larger  angle  as  it  leaves  the  nozzle.  As  a  result,  the  strength  of  the 
impingement  shock  increases.  Once  the  impingement  shock  is  sufficiently 
strong,  the  resultant  adverse  pressure  gradient  causes  a  fraction  of  the 
exhaust  flow  to  be  turned  upstream  (blow-by  flow). 

Theoretical  and  experimental  values  of  the  secondary  flow  rates 
are  presented  as  a  function  of  the  stagnation  pressure  in  Figs.  13  and  14. 

For  all  three  configurations,  theory  provides  a  good  prediction  of  the 
stagnation  pressure  for  which  there  is  no  flow  in  the  annular  gap.  How¬ 
ever,  the  theory  provides  reasonable  estimates  of  the  secondary  flow  rates 
only  for  relatively  low  flow  rates.  This  is  evidenced  by  the  comparison 

between  data  and  theory  for  the  C3/L3  configuration,  which  deteriorates 

6  *? 

as  p^  increases  above  6  x  10  N/m'. 

This  result  should  not  be  surprising.  Furthermore,  as  has  been 
often  noted,  the  theoretical  model  of  the  shear-layer  does  not  reflect 
the  interaction  between  the  approaching  free  shear  layer  and  the  pressure 
rise  of  the  impingement  shock-wave.  Both  of  these  approximations  con¬ 
tribute  to  the  limitation  of  the  flow  model  to  the  prediction  of  the 
lower  secondary  flow  rates.  As  the  model  is  presently  formulated,  the 
secondary  flows  are  only  carried  between  streamlines  within  the  linear 
mixing  profile.  Thus,  the  model  is  not  applicable  to  high  secondary  flow- 
rates. 

For  the  C5/L1  configuration,  the  theory  provides  a  reasonable 
approximation  of  the  crossover  value  of  the  stagnation  pressure  but  does 
not  come  close  to  the  blow-by  flow-rates,  as  shown  in  Fig.  14.  Part  of 
the  difference  may  be  experimental.  Since  the  diameter  of  the  pitot 

probe  was  approximately  50  of  the  width  of  the  annular  gap  for  the  LI  i 

tube  configuration,  inaccuracies  in  the  total  pressure  measured  by  the 
pitot  probes  may  contribute  to  significant  experimental  errors.  Neverthe¬ 
less,  the  theoretical  flow  model  provides  a  reasonable  estimate  of  the 
stagnation  pressure  required  for  zero  secondary  flow. 
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Concluding  Remarks 


An  engineering  model  has  been  developed  to  describe  the  impinge¬ 
ment  flow  field  which  is  produced  when  an  underexpanded  supersonic  nozzle 
exhausts  into  a  constant-area  tube.  Solutions  of  the  theoretical  flow- 
fields  have  been  computed  from  the  transonic  region  of  the  nozzle  throat, 
through  the  nozzle,  and  through  the  impingement  of  the  exhaust  plume  on 
the  launcher  wall.  The  theoretical  solutions  have  been  compared  with 
experimental  data.  Based  on  the  correlations  between  the  data  and  the 
theory,  the  following  conclusions  are  made. 

1.  The  correlations  between  the  theoretical  and  the  experimentally- 
determined  flow-fields  in  the  nozzle  and  of  the  free  exhaust  plume 
indicate  that  the  nozzle-exhaust  flow-field  can  be  well  modeled 
analytical ly. 

2.  With  the  annular  gap  sealed  by  an  obturator  ring,  the  theoretical 
model  closely  predicts  the  static  pressure  in  the  "dead-air"  base 
region.  The  theoretical  and  the  experimental  static  wall-pressure 
distributions  in  the  impingement  region  are  in  reasonable  agreement, 
with  the  best  correlation  for  the  relatively  large-radius  tube. 

3.  For  the  vented  annular-gap  configurations ,  the  theoretical  wall-static 
pressure  distributions  in  the  impingement  region  again  are  in  reason¬ 
able  agreement  with  the  data,  providing  the  impingement  shock-wave 
can  be  calculated  using  the  weak,  oblique,  shock-wave  relations. 

4.  The  theoretical,  flow  model  provides  a  reasonable  estimate  of  the 
reservoir  (stagnation)  pressure  for  which  there  is  no  secondary  flow 
in  the  annular  gap.  However,  the  theoretical  flow  model  represents 
the  actual  flow-field  only  when  the  secondary  flow  can  be  accommodated 
within  the  linear  mixing  profile. 
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(b)  System  produces  nominal  blow-by  flow 


Figure  1.  -  Sketch  of  launch  tube  flow-fields  that  produce  relative 
small  secondary-flow-rates. 
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Figure  2.  -  Sketch  of  impingement  fl 
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Figure  2.  -  Concluded. 
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Experimental  values 
Theoretical  values 


Figure  5.  -  Sketch  of  the  static  pressure  distribution  for  the  divergent 
of  the  C5  nozzle. 
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-  Radial  pitot-pressure  distributions  for  the  exhaust  of 
the  C3  nozzle  (p^  =  8.715  x  10^  N/m?). 
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Figure  6.  -  Concluded. 


Experimentally-  determined  plume  boundary 
Theoretical  plume  boundary  ,  Method-of  -  Characteristics 
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Fiqure  7.  -  Plume  boundaries  for  the  C3  and  for  the  C5  nozzles  exhausting  into  quiescent 
air  (p.,  =  8.715  x  106  N/tn2). 
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Figure  8.  -  Base  region  pressure  for  a  sealed  annular  gap. 
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-  Static  wall-pressure-distributions  for  C3  nozzle 
exhausting  into  a  launch  tube  with  a  vented  annular  gap 

(8.715  x  106  N/m2). 
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